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ABSTRACT 


A photometric re-examination of the eclipsing variable 
EI CEPHEI was undertaken using the 12.5 inch Cassegrain 
telescope housed at the Devon observatory. The intensity 
measurements were made in four colours defined by the uvby 
filter system. The method of analysis employed to solve for 
the -vacity, “luminosities, and angle of inclination of the 
orbit was the procedure suggested by M. Kitamura in 1965. A 
FORTRAN program, EVAL, was written to calculate the various 
parameters required to use the Kitamura method. Observations 
of the variable were made from August, 1972 to October, 
1973 from which 208 data points in each colour were obtained. 
The data when plotted to form the light curves demonstrated 
a high degree of scatter. Although it was possible to 
easily observe the minima, the light curves were of insuf 
ficient acciiracy to yeild a final solution for the orbital 
elements. It was, however, possible to show that the 
Cramary Minimum we. cue, to.an occultation (eclipse or the 
smaller star by the larger component) while the secondary 


minimum is due to a transit (eclipse of the larger star). 
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CHAPTER. L 


INTRODUCTION 


Eel). BINARY STARS 

With the development of the telescope it was found 
that certain stars formerly observed to be singular could 
be resolved into components. Such stars were considered as 
little more than astronomical curiosities since it was pre- 
sumed they were merely stars that possessed similar spatial 
coordinates. The attitude towards these objects changed 
however with tne announcement by Halley in 1718 that some 
of the brighter stars, Sirius, Arcturus, and Aldebaran were 
in motion. In 1767, John Michell proposed an argument 
based on the theory of probabilities that some of these 
double stars could be systems joined together under some 
general law. The real beginning to double star astronomy 
however is credited to Sir William Herschel, who carefully 
investigated such objects throughout the last quarter of 
the eighteenth century. The culmination of his work was 


(2) 


presented in a paper written 1803 in which he proved by 
actual demonstration from measurements he had made that 
certain double stars were true binary systems. Stars 


where it is possible to optically resolve the individual 


components are designated as visual binaries. 
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Another class of binaries was discovered with the 
aid of the spectroscope. Stars that still appeared single 
in a telescope were found to display periodic oscillations 
of spectral lines about their mean positions. This movement 
of the spectral lines of the so-called spectroscopic binar- 
ies was attributed to the Doppler effect resulting from the 
motions of the stars in their mutual orbits. If the plane 
of the orbit of one of these spectroscopic binaries should 
happen to fall close to the line of sight from the earth 
then it is possible that the two components will periodi- 
cally eclipse one another. During the passage of one star 
in front of its companion, the apparent brightness of the 
system as a whole will diminish. Such a celestial object 
is classified as the eclipsing binary or eclipsing variable. 
A common point among binaries is that by proper mathemati- 
cal analysis, pertinent information about the star system 
can be obtained. This is especially true for the eclipsing 
binary and therefore it has been chosen as the basis of 


this thesis. 


Leas PHOTOMETRY: 
The distinguishing feature of eclipsing binaries is 
their light variations. Therefore, it is essential as a 


method of detection to measure any fluctuations in light 


do 
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intensity. This problem can be resolved by employing the 


general concepts of photometry; astronomical photometry is 


concerned with measurement of stellar brightness (magni- 


tudes). The first device used as a photometer was the eye. 


Stars were compared to cnosen standards in an attempt to ca- 


talogue them according to their apparent brightness. This 


method proved extremely inaccurate due to the eye's inabil- 


Ley 
eye 
was 
the 


the 


to accurately evaluate intensity intervals. To aid the 
in making more accurate estimates the use of null-methods 
instituted. One such device applying this technique was 


(2) 


polarizing photometer. Here a telescopic image of 


real star was compared to the image of an artificial 


star whose brightness could be varied precisely by a rotating 


polarizing prism. Although it was possible to obtain stellar 


magnitudes with uncertainties as low as 0.04 magnitude, there 


still existed large systematic errors due to differences in 


color sensitivity between different observers. 


With the development of the photographic plate in 


1850 a new system of measurement became possible. The fun- 


damental principle of photographic photometry is that equal 


intensities produce equal photographic effects under identi- 


cal 


conditions. 
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These conditions included the emulsion, the exposure time, 
temperature and humidity during exposure, and the size and 
structure of the optical image. Because of the difficulties 
in controlling changes in some of these conditions and the 
non-linear response of the emulsion with exposure time, 
photographs of the brightest and faintest stars are often 
incorrectly exposed. Despite these drawbacks, photographic 
photometry has the unique capacity to simultaneously record 
large numbers of objects in a limited field which is 
indispensible in stellar statistics, variable star surveys 
and investigations of rich star clusters. By establishing 
reliable magnitude sequences in the region under study it 
became possible to reduce the errors in estimation to 0.02 


magnitude. 


The most recent advance, and that used Peonehoue 
this project, has been the development of the photoelectric 
photometer an electronic device that displays a linear 
response over a broad range of stellar magnitude. The basis 
for such an instrument is the photoelectric effect dis 
covered by Hallwachs in 1888. Hallwachs observed that when 
light fell on certain substances, a measurable electric 
Current proportional, to the incident light intensity resulted 


and, therefore, an accurate measure of this current could 
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produce an accurate determination of the light intensity. 
Early work in this field of photometry was carried out 
independentlyiby JieStebbinseinethesUvserand!P a Guthnik@in 
Germany. The first photocells used had a selenium photo- 
cathode of low sensitivity which restricted work to the 
very bright stars. After 1911, more sensitive potassium 
hydride cells were constructed. Despite the development of 
these more sensitive cells, the measurement of the output 
current still required the use of elaborate electrometers 
which consumed valuable quantities of man-hours in the 
course of measuring a single light-curve. With the improve- 
ment of electronic amplifiers, it became possible to in- 
crease feeble cell currents to the point where ordinary 
galvanometers could be used. A major improvement occurred 
with the production) of, ‘the: RCA» 1P21° photoelectric cell’ en-— 
ploying the use of secondary emission which coupled a high 
internal gain with a broad linear response. This new cell 
together with new recording instruments resulted in accur- 
acies of 0.01 magnitudes -<Theslimiting»effectsvon the ac= 
curacy of photoelectric measurements are imposed by the 
atmosphere (e.g. changes in transparency and scintillation 
effects); and noise in the photomultiplier) which) will lin= 


it the magnitudes that can be measured. 
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1.3. THE CLASSIFICATIONS AND ASTRONOMICAL IMPORTANCE OF 
ECLIPSING VARIABLES 

By using the photoelectric photometer and its an- 
cillary equipment, it is possible to obtain highly accurate 
observational data. These data for an eclipsing binary con- 
Sists of a relation between apparent brightness and time. 
A graph of these measured quantities is called a light- 
curve. A study of the many light-curves obtained to date 
has shown that it is possible to classify the eclipsing 
binaries as follows: 

(1) no evidence of any significant extra-eclipse 
variation, therefore the light between eclipses 
can be regarded as constant. 

(2) theriight isefotindstoevany duc.to, retlection 
and ellipticity effects but the moments at 
which the eclipse begins and ends are still dis- 
tinguishable. 

(3) the changes in light between eclipses are com- 
parable to those caused by the eclipses and 
make it impossible to approximate the beginning 


and end of an eclipse. 


These three classifications which are based on the 
degree of proximity effects correspond approximately to the 


following models which are named after the prototype star: 
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(1) Algol: component stars well-separated, so 
Stars are essentially spherical; orbit 
is very close to circular. 

(2) 6 Lyrae: components so close as to be vir- 
LUaliv in contacts, tidal efrects ai6— 
tort stars into elongated ellipsoids, 
eccentric orbits. 

(3) W Ursae Majoris: the components are physically 
united and as such act as a singular 
unit, period of revolution is less than 
12 hours due to compactness of system. 

These classifications, corresponding to theoretical 

models are important in the analysis of all light curves. 
It is possible through techniques to be discussed in the 
next chapter to solve for: 

(1) the radii of the two components rye ro in terms 
of mean distance, a= (a, + ay) 

(2) the luminosities Lys L. in terms of the total 
luminosity of the system outside eclipse. 

(3)7 the angle of inclination Of the orbit, 1eco the 
line of sight of the observer. Angle i ap- 
proaches 90° as the line of sight approaches 
plane of orbit. 

(4) the ratio of the surface brightness 
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ifathepeclipsing binaryeis:also studied spectroscopically 
and if both spectra are visible, then an analysis of the 
velocity curve given by the Doppler shift of the spectral 
lines will yield the quantities ay Santi, as sin i and my 
ke. m5 eine te Therefore a combination of photometric 
and spectroscopicedata will yield» the absolute: values of 
the radii, masses, and densities of the eclipsing system. 
When the problem of an eclipsing binary system is studied 
astrophysically more information than merely the dimensions 
or the surface characteristics of the components or the 
geometry of the orbits are provided. With a single star it 
is impossible to probe the gravitational field or learn any 
thing about the surface brightness. However, place a com- 
panion in close proximity and properties of. the combined 
gravitational field can be deduced from characteristics of 
the motion (e.g. apsidal advance) or by their mutual dis- 
LOLreLon resulzingyan chengésiinetheslight tusverii The diss 
tribution of intensity over! the apparent) disks! canbe 
studied through®thepvardation an brightness provoked by ax— 
jal rotation of the distorted components. Non-radial 
Gecililations of. tidal origin.in.a close. eccentric system 
can (through tidal lag) disclose the degree of viscosity 


of the outer layers of the component stars. 
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1.4 STATEMENT OF THE PROBLEM 

The basic aims of this research project are listed 

below: 

(1) Obtain complete photoelectric light curves in 
the uvby intermediate-band photometric system 
following the procedures of differential 
photometry. 

(2) If possible, reduce the photometry to the 
standard uvby system. 

(3) sin carrying out (1) and (2) it willebemeces- 
Sary to determine extinction coefficients. 

(4), “In ‘carrying out. {1) and) (2) 16 wiil ber possibie 
to formulate conclusions regarding the preci- 
sion with wnich intermediate-band photometry 
can be done on faint stars with a small tele- 
scope and an unrefrigerated photomultiplier. 

(5) eAnalyse the light curves for the orbital ele= 
ments and combine these with published spec- 


troscopic elements. 


The eclipsing system under investigation is EI CEPHEI 
[BDi +h 759079 =eHD205234.5+6C30132)..,iThisystar,is,a Henry 
Draper spectral type FO with a visual magnitude of 7.67. It 


is of particular importance since EI CEP is one of only 19 
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known main sequence eclipsing variables of spectral classes 
B9-G for which accurate determination of the component mas- 
ses Should be possible. EI CEP was first established not 

as an eclipsing variable but as a double-line spectroscopic 
binary by J.S. Plaskett of the Dominion Astrophysical Obser- 


(3) 


vetonw(YLetoria) an 971, A provisional set of elements 
Of the spectroscopic orbit was calculated in 1935 by W.8. 
Harper of the DeATOu The possibility that EE CEP was 
also a photometric binary was suggested by WwW. Strohmeier ‘°) 
at the Bamburg Observatory in 1958 during a search for new 
variables. A photometric study was undertaken by A. Abrami 
of the Astronomical Observatory of Trieste during the per- 
iod of 1959-1965 using an unfiltered photometer. Abrami's 
results confirmea Strohmeier's earlier suggestion that EI 
CEPHEI was indeed an eclipsing binary and he further cal- 
culated it had a period of 8.439365 days with two narrow 
Minima of approximately the same duration and with depths 
of about 0.507 and 0.422 magnitude respectively. The anal- 
ysis of Abrami's data for the orbital elements was hampered 
by the presence of third light constituting 9.0% of the 
total luminosity. The inclusion of this extraneous light 


brought with it a measurable uncertainty in the values of 


the photometric elements. In 1969 S. Kriz ‘/) of the 
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Czechoslovak Academy of Sciences used the combined photo 
metric and spectroscopic data for a comparison with the- 
oretical models of main-sequence stars as proposed by Iben 
and Horn. Kriz concluded that the models agreed well with 
most of the observational data available except for EI CEP- 
HEI which deviated in all comparisons. This deviation from 
the theoretical predictions was attributed by Kriz to the 
considerable uncertainty arising from the combination of 
Harper's preliminary spectroscopic data with Abrami's some- 
what uncertain light elements. Kriz finally concluded that 
these uncertainties were large enough to warrant a new set 
of orbital elements. The first reinvestigation of EI CEP 


8) 


waS carried out by D. Popper ‘ in 1970 with a recalcula- 
tion of the syectroscopic orbit from a new set of spectro- 
grams. The objective of this thesis will therefore involve 


a multi-colour photometric study of EI CEPHEI in an attempt 


to obtain more accurate photometric elements. 
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Harper (1935) 
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SPECTROSCOPIC ELEMENTS 


Popper (1970) 


Period (days) 84455 8.439365 
Eccentricity 0.039 

K, (km/sec) (nel Ble 0s 
Ky (km/sec) 80.28 76.9 ¥ 065 
To 2417978327 2436820 .4665 
Bisinwa(solar radii) © 26.7 ZOdiaa 1 
m1 sin i (solar masses) 1.80 incGiaa 2003 
n, Le i (solar masses) 1.79 ineycie 2003 


PHOTOMETRIC ELEMENTS 


Abrami (1966) 


k (vr, + rg) 0.774 
Inclination, i (degrees) 87.97 
r Os067 
Ss 
ve 9 Pe a 
& 
fie 0.385 
Ly Os Were 
0.093 
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a, (kn) ga2s 10” 
a, (xn) 9.32 x 10° 
m4 (solar masses) dea AS 
n., (solar masses) HLT AS) 
x (solar radii) PES 
lp (solar radii) a0 


The quantities Ky and Ky represent the semi-amplitudes of 
the radial velocity changes versus the phase (position of 
the two components) for Star 1 and Star 2. ‘The value 
denoted by the letter a is the apparent separation of the 
two stars expressed in solar radii. To represents an 
accurately measured time of primary minimum, The values 
Peilgyly=ls,lp=h, represent the radius and luminosity of 
the small and large component respectively normalized by 
a (=a, +25) for the radii and the total luminosity outside 


eclipse for Ls and L The angle of inclination is the 


g&° 
angle formed by the line of sight of the observer and the 
normal to the plane of orbit of the system. The presence 
of third light indicates that an unseen object is contri- 


buting to the total output of the system. 
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CHAPTER 2 


THEORETICAL TREATMENT OF ECLIPSING VARIABLES 


200) GEOMETRY OF THE ECLIPSES 


The interpretation of any eclipsing star system is 
based on a comparison of observations with a simple model or 
a series of increasingly more complex models. In the case 
of an eclipsing binary, the simplest system consists of two 
spherical stars of uniform brightness moving in circular 
Orbits. Although this “model is not entirely realistic it 


does provide a useful and often accurate starting point. 


At the conjunctions of the stars there will be alter- 
nating eclipses, while out of eclipse the light output will 
remain constant. Because of the spherical nature of the stars 
and their circular orbits, the minima will be symmetrical, of 
equal duration but (generally) unequal depth, and will occur 
exactly one-half an orbital period apart. A careful study of 
the light curve will reveal either a continuous variation of 
light during both minima corresponding to partial eclipses or a 
phase of constant light at the bottom of one (or both) minimum 
due to alternating total and annular eclipses. In general any 
Minimum, total Or partial, arising from an eclipse of Che smal- 
ler star by the larger component is termed an occultation where- 


as the reverse configuration is called a transit. Difficulties in 
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discriminating between the various types may arise if the 
eclipse ts nearly grazing or if only one minimum has been de 
tected. In studying a light curve where both eclipses have 
been observed, the deeper minimum is referred to as the primary 
and the shallower minimum as the secondary. The deeper min 
imum will always correspond to the eclipse of the star of 


Greater suriace brightness, 


At any moment during an eclipse, one component eclipses 
a certain area of the other; however, during the corresponding 
phase of the alternate eclipse, the geometrical relation of 
the two disks remains unchanged and the other star eclipses an 
equal area (FIG. 2.1) but not necessarily an equal fraction of 
the total area of the other star. If, as we have assumed, 
both components are uniformly bright, i.e. the distribution of 
surface brightness is equal at each point on the stellar disk, 
then the changes in light during both minima should be similar 
and differ only on a scale corresponding to the ratio of the 
surface brightness of the two stars. This type of simple 
relation is complicated by the fact that most stars display 
some degree of limb darkening. This phenomenon was first noticed 
from observations of the sun. The limb (thee-Outer-most part 
of the visible surface) of the sun was observed to be darker 


than the central region. The reasoning behind this effect is 
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GEOMETRY OF THE ECLIPSES 


The eclipse is an occultation if the eclipsed (shaded) area 


is on the smaller star. 


The eclipse is a transit if the eclipsed area is on the 


larger star. 


Y 


a 

oy st 
a -_: 

" | Tete oselfae at? 26 ar" 

ae aa 


. | ; 
gwen (Webra) -bedgtion ads Wi coteetiaids -& of ange 


wl? J Gi oe » Seeqlizue. § 7 s 24 panz a af cagifes 
ee _ _ a 


tera teygel 
o d 


Bo) 


These conditions included the emulsion, the exposure time, 
temperature and humidity during exposure, and the size and 
structure of the optical image. Because of the difficulties 
in controlling changes in some of these conditions and the 
non-linear response of the emulsion with exposure time, 
photographs or the brightest: -and®faintest®starsrarevorten 
incorrectly exposed. Despite these drawbacks, photographic 
photometry has the unique capacity to simultaneously record 
large numbérs of objects in a limited field which is 
indispensible in stellar statistics, variable star surveys 
and investigations of rich star clusters. By establishing 
reliable magnitude sequences in the region under study it 
became possible to reduce the errors in estimation to 0.02 


magnitude. 


The most recent advance, and that used throughout 
this project, has been the development of the photoelectric 
photometer an electronic device that displays a linear 
response over a broad range of stellar magnitude. The basis 
for such an instrument is the photoelectric effect dis 
covered by Hallwachs in 1888. Hallwachs observed that when 
light fell on certain substances, a measurable electric 
current proportional to the incident light intensity resulted 


and, therefore, an accurate measure of this current could 
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associated with the varying optical depths of the sun's outer 
surface. If we observe the sun near the centre we receive 
light from a layer corresponding to a certain depth and temp- 
erature. However, near the edge one can see a comparable 
distance but the deepest layer is higher than at the centre 
(seemEIGt 222)3 therefore; othe temperaturenis owerfandethe 

sunt appears darkeroat that point.) SThe net effectrofelimb-dark— 
ening» on the light curve is to decrease the overall loss of 


light during an eclipse. 
PLAS THEORETICAL LIGHT CHANGES 


For our presumed model involving circular motion, the 
longitude of the orbit, 8, measured from primary minimum is 
given by 


6 = 360~ (t-t) 
P 


where P is the period (time between two successive primary 
minima) and te the epoch of primary conjunction both of which 
are expressed in Julian days. In order to relate the position 
as a function of time to the light change during an eclipse, 
we consider two stars with an apparent separation 6 between 


their centres where from geometrical considerations: 


— = egos! + ein“isin-s (201) 
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SCHEMATIC EXPLANATION OF LIMB DARKENING 


The light from the centre is emitted from a deeper region 
than at the limb. The angle formed by the intersection 
of the surface normal and the line of sight is referred 


to as the angle of foreshortening. 
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The fractional loss of light at any moment during an eclipse 
can be defined as a function depending upon ra the radius 
of the smaller component, Dye the, radians ofethedlargergstar, 
and,6o, theyeseparation ofether centres. -\Forethescasei ofeuni- 
form, brightness, sthesfractionaladossief lightyaf,iasetheira- 
tio of the eclipsed area to the total area of the eclipsed 
star, while for the completely limb-darkened case, f is the 
ratio of the light emitted from the eclipsed area to the to- 
tal “bight of the star being eclipsed. The function £ is 


homogeneous in the three elements ror x andaoy I tkcan 


jay! 
therefore be made to depend on only two ratios formed between 
them. The most advantageous of these are the ratio of the 


radii, k and the geometrical depth of eclipse, p, where by 


definition: 


sets ne O<k<l1l (222) 
Ws 
and 
p= aa ; p = 1 monent of first contact 
rs p = -l moment of internal tangency 


The two radii are normalized to the unknown radius of the 
circular orbit. The loss of light: during eclipse can be 


expressed as (9). 
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L = | COS ey -oG (223) 
S 


where J is the surface brightness, y the angle of foreshort- 
ening, do the surface element, and s the region of integra- 
tion. From the principles of radiative transfer in stellar 
atmospheres, the change in surface brightness across the disk 


Or a Star follows a cosine Jaw of the form 
J(Y¥) = J(o) (l-utu cos y) (Zea 


where u is the coefficient of limb darkening and y as before 
the angle of foreshortening. If the change in light as 
found from equation 2.3 is compared to the light emanating 
from the entire disk of the eclinsed star, the resulting 
value will represent the fractional loss of light where the 
f-functions in general are complex mathematical expressions 
involving transcendental functions. The fractional losses 
of light have been evaluated for the uniformly bright star 
(u=0), £ and the completely limb-darkened star (u=1), c?. 
such that for arbitrary limb-darkening, the loss of light 


can be expressed as a weighted mean of the two extreme situa- 


tionse(o)}:: 


fea oii yer He 2ue E (225) 


nN 


[4 


eee ey Saved) «a - 


» : 7 
<#2exies qo te et pite oa? 4 saoareriioizd sontuls sly of Us 
~auyedat 43 mises oe aie .saetel> @odtave si Hy 
aeiiase 64 tatengy? evdantees do molqyioniny, o8F was 7. 


s 


weit) hf Gengee ceosifei eh wiatorc cc aumadl> ahs | 
<7101 G3 8m wa! sateor + etafias 1980 BOE 


; . 
(Red! (vy; 260 21-2) (a) © gerh : 


acing as 7 lau akeugiet.' dald ic sate steoo ene abe etity: 7 
. eA deli AS meneD At? TT). Holes tutteertd “20 apap Oe 
end se; 20 ones gia at Waudtrd «<c '.° @otsse7e we eae o 

yinastuent Of, Wode Papglian ast Te ABels wlohe ed3 a 
eas wzaciv. ciithd Yc audl Ladedsonx? an2 Jao te lyee Oh Gale 
enviversias iethtmeeiinn +4i\yro bat, toterny of. eeehson0t-D 
eace5! Léadi2oe1? | of réot!amt™ Lasovtunscanedspabvioval 
gaze Shp) sd yirrest din Ne G3 “tshiluve ase evar siest — ; 
Ys » he) ante benvAteesdatl yipfulemoo srt tae, ty ere 
Spit Yo sco] ai \patanstededuti yasitidws xh Gade dots 
wepdle anecsxe Ws. oF> ID mage Latipias « we bommumsque sel neo 
rte} bs 
(a2) Sg a. My tu ane = 3 
rE 


Thereforesit isepossiblento calculate the: fractionalylosses 
of dighteyfrom tabillesiot £", ? which have been extensively 
computed by Zessevitch (10) and Merrill (11). If we consider 
the individual luminosities of the two components normalized 
by La + Ly = lJ, then the foss of light during any minimum can 


be computed from (9): 


a ~ 
Se ie) as 
= 1-f (kK p)L. 
i=l 
oe dee Tage ii (2-6) 


depending upon whether the eclipse is an occultation 


Om .actransit. 


2.5 tHE ANALYSIS OF A LEGHT CURVE 
The’ most diificult and basic problem confronting the 
analyst is to separate the light changes due to proximity 


effects from those due to the eclipses. However, in order to 


2 


treat the method of analysis in an abbreviated form, consideration 


will be given only to systems of spherical stars sufficiently 


well separated as to introduce no auxilliary effects. 


The fundamental problem of the eclipsing system can be 


reduced to two basic equations that under appropriate analysis 


will yield the unknown elements. 
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The first such equation arises from equating the value of 
Gsasegiven by “equation 2.1 with that given in equation 2.2. 


The result being 
ry (1+kp)? is aeee sin? a cos*i (2a 


The second equation is found by comparing the normalized 
light loss (the fractional losses of light £U, ae normal- 
ized at the moment of internal tangency, p = -1) at differ- 


ent phases with the value of the depth of minima such that 


Ok, p) = H(k,p) = 1-2 aU P)) = no 
E(k 1) 1-A 


where n = 1-2 
- 


ee 


and i} is the luminosity at midminimum. The quantity 

a = pikrD @ is the maximum obscuration i.e. the maximum 
value of the normalized light loss at mid-eclipse while a. 
is the value of the geometrical depth of the minimum also 
at mid-eclipse. 

THe L£irst. attempt @t a direct solution for the orbital 
elements was made by E.C. Pickering (12) on Algol (8 Persei). 
Despite Pickering's efforts, the problem remained a chal- 
lenge until Henry Norris Russell proposed (13) an efficient 


and workable method in 1912. 
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@) RUSSELL'S METHOD 
The first fundamental equation 2.6 in the case of a 


total eclipse contains three unknowns k, r,_, and i which 


b 
should be soluble by choosing three points on the light 


curve. Russell's idea was to eliminate ry and i, then 


solve for k as the root of the determinental equation (9): 


ee ; 2 

Sin” (1+kp,) i 
ee 2 

sin” , (1+kp.) iL = 0 
ee 2 

Sin” 3 (1+kp 3) a 


Russell, in order to solve for k, assigned two of the points 


as specific values (a, =" 0.6, a5 = 0.9) while retaining the 


third as arbitrary. Expanding the determinant reveals that 


sin?@ = A + BY(k,p) (2.8) 
A = sin“@, 
B = sin76, ve eames 
re 
Bpep we kip ps) 
p(k, o) Fe ae UNNe 
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The whole nature of the problem was then confined to the 

function ~(k,a), but this was resolved by tabulating values 
for various degrees of limb darkening. Once the value of k 
had been chosen so equation 2.7 had been satisfied, the re- 


maining elements were provided by: 
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where >, (kK), 5 (x) were again tabulated functions which 
facilitated a more workable solution. The concept of tabi 


lated functions originated with the Russell method and be- 
cause of its usefulness was used by most subsequent investi- 


gators. 


Despite the overall success achieved by Russell's 
method, certain defects existed centered around the use of 
the pivotal points... By selecting specific points, arbitrary 
weights were assigned to different parts of the light curve 
which made it impossible to ascertain the errors in the 
solution. To overcome this problem the concept of the 


iterative method was developed by Zdenek Kopal (14). 


(ite) KOPAL'S METHOD 

Again as a starting point, Kopal chose equation 2.6. 
At the moment of inner contact p = -l, the above equation 
reduces to 


sine. sini + cos*i =r ete Key (2.9) 
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By SOlving for sine 8 from equations 2.6 and 2.8, Kopal 


obtained the equation (14) 


(p*-1)c, + DUS hee ees retina (UE 
oy ty 2% 
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Using a preliminary value of k and tables of p(k,a) equation 
2.9 could be solved by least squares using each point on the 
light curve. A comparison of the initial value of k to k = 
C,/C. determined whether or not another iteration need be 


performed. If not one could solve for the unknowns from 


r= (Cjsin i 


ry es C5sin al 
NiCy 
2 (C1 2eS) Fae 
cot i= Zee 3 
Sy 


The use of Kopal'’s technique as a by-product of the least 
squares analysis provided the weighting of each element and 


the probable errors incurred. 
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2.4 METHOD OF KITAMURA 


The preceding procedures for the analysis of light 


curves proved disadvantageous in two areas. These were: 


(1) it was impossible to know prior to the evalu- 
ation* o£ the-elements-which of the»two minima 
was due to an occultation and which was due to 


ar transi. ts 


(2) the method of analysis was different for the 
different types of eclipses (partial, total, 
annular) which in themselves could be difficult 
to distinguish by mere inspection of the light 
curve. 

These difficulties were overcome by the recent development 

of a new method as suggested by Masatoshi Kitamura (15,16). 

This method of analysis is dependent upon certain characteristic 
functions of the eclipse derived from integral transforms 


of the ascending and descending branches of each minimum. 


As in previous methods, the distribution of surface 
brightness is assumed to follow the cosine law given by 
Smeets ETheduminosatyootethessystem during an eclipse is 
given by 


&(0) = 1-L,£(8) eae) 
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where as before the fractional loss of light is given by 
equation 2.5. Two types of integral transforms, the incom 
plete Fourier transforms are defined as (15) 

ak 
S_= of (8) sinnédé 


n 


XX 


hore 
‘ee -f ok (6) cosné6dée 


@ is the phase angle of an observed point outside of 
eclipse and n is zero or an integer. Another pair of trans - 


forms can be Lormulated as (25) 


6 
s =I £(9) sinn6édg 


Cc ace £(6) cosnéede 


es is the phase angle at the moment of first contact given 


by 


serail = ses -1l 
The. fetjo 2 2 
tL (ro + ty ) 


By substitution of equation 2.11, a relation gan be found 
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Lis, = i089 - = Ae et Oe ES ie eee 
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The conversion of the observed light curve into its incom 
plete Fourier transforms can be performed numerically with 
the aid of an electronic computer. The problem of finding 
a solution reduces to the simultaneous determination of the 
geometrical elements of the eclipse rar XH and i. In order 


to solve for these three unknowns, three independent relations 


are required (15): 


(T)° ratio between’ i. sey, > Loco" for speci fied 
art 1vn 
values of n 
Fn i. len - na 
riilas <n 
(2) the ratio between Lic, and the depth of 


Minimum, A, 


(3) combination of both minima 
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Ds = (Ljc,)p p = primary 
DA (Lie )s S = secondary 
solving for luminosities 
Bo =_D 
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while on (Cede are given by theory. 


The three functions Poe Ey om are now considered as 
the characteristic functions of the eclipse. The solution 
oO oT the ang 191n) principle Loviows, from es compari sonror 
the values of the functions E and oe deduced from the observa - 
tions for either minimum with their theoretical values as 
functions of the geometrical elements and use of the third 
relation to establish the best set of elements. The same pro- 


cedure is applied to any type of eclipse. 
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The type of eclipse occurring at each minimum is obtained 


byecons i deration Of the theoretical Natureso. the eclipse 


and the mathematical properties of the functions discussed 
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in section 2-2. Kitamura has shown that 
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where. the equal sign holds-for 
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(1) grazing eclipse 


(25) u_ = Uy, = 0 


Because of the advantages enjoyed by this method and its 
application to computer analysis, it was chosen as the 


procedure to be used in the examination of EI CEPHEI. 
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CHAPTER 3 


OBSERVATIONAL EQUIPMENT AND EXPERIMENTAL TECHNIQUE 


3.1 OBSERVING EQUIPMENT 


The observations of EI CEPHEI were carried out at the 
Devon Observatory located approximately twelve miles from 
the city of Edmonton, Alberta. The equipment housed there 
and used throughout the study of EI CEPHEI includes: 

I2% reflecting telescope and single channel photometer. 

high voltage d.c. voltage power supply 

digital voltmeter and didqitel ‘clock 

aucLomaticerecorder —- printer 

multi-band radio 
The accuse ltion oF the ~lour=band sadio in early 1973 per 
mitted monitoring of the digital display clock by comparison 
with transmitted radio timing signals. The physical construc 
tion of the photometer was altered slightly during the course 
of observing. It was noted that different readings were 
being obtained when the telescope was shifted to the opposite 
side of its mounting pier. This effect was attributed to a 
sagging of the photometer in the collar used to couple it 
to the telescope. To rectify this situation, three supporting 
members constructed from metal rods were attached to the 


photometer and telescope to ensure a rigid structure. 
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(ze) TELESCOPE 

The telescope used throughout the observing sequence 
was a Cassegrainian - type reflector equipped with an equa- 
torial mount. This instrument was obtained in 1967 from 


Tinsley Laboratories, Berkeley, California. 


The primary mirror is 12.5 inches in diameter with a focal 
length of 48 inches while the secondary mirror is 3.75 inches 
in diameter. The initial installation of the equipment 
brought with it certain mechanical problems. A close examin 
ation of these problems and methods used to correct them are 


described by Bo Brown (17). 


(11). PHOTOMETER 
The basic components constituting the photometer are: 
field eyepiece and assembly 
diaphragm 
diaphragm - viewing telescope and assembly 
filter assembly 
Fabry lens 


photomultiplier 


An outline of the basic construction and a description of the 
individual components and their use is again given by B. Brown 


(37) F@eAlthnough the photometer 1s "essentially the same as” that 
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described in the above reference, certain changes have 


been made. 


(a) photomultiplier 


In previous observations with the photometer an RCA 
or Sylvannia 1P21 photomultiplier tube was used. This cell 
was characterized by an entrance window on the side which 
meant the photomultiplier had to be mounted perpendicular 
to the main photometer housing. In its place was substituted 
an E.M.I. 9502S end-on photomultiplier cell which can be 


placed directly into the photometer housing. 


Cis) FILTERS 
Previous multicolour photometric work at Devon was 
carried out in the UBV filter system developed by Johnson 


and Mergen oe 


Observations on EI CEPHEI were done, however, using the 
intermediate band filters of the "uvby" system. This system was 


(18) who’ describes Tt as "an extension 


developed by B. Stromgren 
to the MK (Morgan-Keenan) classifications and UBV photometry. 
The use of intermediate band filters allows a two dimensional 
classification of stars founded upon particular spectral 


indices which are based on the measures of selected spectral 


features. 
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Intermediate band photometry is characterized by filters 


fe) O 
WowCcnenave-oahalt—widtheO ee 00 AL 300 Ag sinsparti cular, 


the filters defining the system in use are: 


(a). u 


(by v. 


(cays oo 


cb a’, 


glass filter (8mm UG11 + 1mm WG3) 
centered about 35004 
transmission region: - above atmospheric cut 
off but below Balmer discontinuity 
interference filter 
centered about 41004 
above region of crowding of hydrogen lines 
as they approach the Balmer discontinuity 
contains the line Hé 
interference filter 
centered about 47000 
contains no strong spectral features 
interference filter 

fe) 


centered about 5500A 


contains no strong spectral features 


Three separate indices have been developed for use with this 


system. They are the colour index (b-y) and the two colour 


qaifferences m4 


(v-b) - (b-y) and aos (G=7)  —s1V—-)).. 
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The colour index (b-y) displays properties similar to the 
B-V of the UBV system but is influenced less by metallic 
absorption lines. The colour difference my is a good indi- 
cator of metal content. The other colour difference c, acts 


ur 


as a good index to the Balmer discontinuity. 


It has been found that the “uvby" system displays 
no effects on band-width such as extinction, transformation 
or interstellar reddening and is also uninfluenced by strong 
spectral features other than HS, Balmer discontinuity and 
the line blanketing. The system is filter defined and hence 
depends little on the atmosphere of the Earth or the photocell 
in use. 


(iv) Miscellaneous 


The remaining equipment was obtained from Hewlett 


Packara and consi cits. Of: 


3340A Digital Voltmeter 
J/4562A, Digital Recorder 


Harrison 6525A DC Power Supply 
(0-4000V, O0-5mA) 


The operational use and limitations of each instrument are 
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. described in detail in the appropriate operating manuals pro- 


vided by the manufacturer. 


B62) OL oLEMATICS 

The first accurate catalogue of stars compiled by 
Hipparchus around 135 B.C. and later refined by Ptolemy to 
form the ALMAGEST listed all objects in five magnitude inter- 
vals. The brightest stars were designated as m = 1 with 
those stars just visible to the naked eye as m= 6. An inves- 
tigation by Sir John Herschel pointed out that the difference 
between a first magnitude star and a sixth magnitude star cor- 
responded to an intensity ratio of approximately one hundred. 
N. Pogson in 1856 proposed a system whereby the intensity ra- 
tio was taken to be exactly one hundred such that each magni- 


tude interval corresponded to (100) °- 


OL 2. 512) in Intensity. 
Pogson's definition together with Fechner's psycho-physical 
law of sensation (sensation <log stimulus) led to the modern 


system of photometry known as the Pogson Scale which states 


m= -2.5 log luminosity me he 


orm, - ™, = De ae OC (2/25). oats) 


A measure of the observed intensities of two stars should 


therefore give an accurate evaluation of their magnitude 


difference. However in reality direct measures obtained in any 
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photometric study are functions of 

(1) parameters related to the equipment 

(2) the atmosphere through which the light must 

travel 

(3) the wavelength of the light! 
The observer must therefore treat his data in such a manner 
as to make allowances for all effects not directly attrib - 
Utable to the object under study. whe reduction of <he data 


consists of two major steps: 


(1) an adjustment for atmospheric extinction 
(2) a transformation to a standard system of 


magnitudes and colours. 


(7) EXTINCTION 

During its passage through the Earth's atmosphere, 
the starlight is reduced in intensity by absorption or through 
scattering processes. The degree of intensity loss will be 
Girecely proportionalsto the thickness of the satmospheres to 
be traversed. The relation that governs the rate at which 


the light is diminshed in magnitude is 
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Here X, the path length, is expressed in units of air mass 
at the zenith of the observer and k, the extinction coeffic ~ 
ient is a measure of the light loss in magnitudes for a star 
also at the zenith. From equation 3.2, a value for the extinc-— 
tion coefficient can be obtained from the slope of a graph of 
observed magnitude versus the air mass. To a first approx - 
imation, a value for the relative air mass is given by the 
secant of the zenith distance, 4%. From geometrical consid- 
erations the value of sec 4% is 

sec % = (sinésingd + cosdcosd cosh) + 
ITP When? io) 1s the observers latitude while 6 and fh are the 
declination and hour angle of the star. If a more exact value is 


ay 


required X may be obtained from the polynomial expansion 


X= "sec 4 — 0.018167 (sec 4 = 1) = 


O;002875 Ssec™a== 15% - 


0200060S3e(sec4 45> ae 


As mentioned before the "uvby" system incorporates three 
indices b-y, M54 and Cy: Each of these values is individu - 
ally affected by the presence of extinction and each must 


accordingly be corrected. The resulting extinction equations 
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From investigations carried out by Crawford and Barnes 
the Kitt Peak National Observatory, it has been found that 
mean values of the extinction coefficients are adequate for 


all but the most uncertain of nights and that no second order 


extinction effects are present. 


(ea) TRANSFORMATION 

The light received by the photometer will be charac 
terized by a certain wavelength and band-width. The response 
of the photometer to this radiation will be dependent upon 
its optical components and the spectral sensitivity of the 


cathode. Therefore each combination of telescope, filters, 


and photomultiplier will define its own natural colour system. 


In order for a set of observations to be of use to another 
obseryer with a different natural system, it is necessary to 
transform all measurements to a universally accepted standard 


system. 
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For ideal black-bodies, the colour indices measured on dif- 
ferent natural systems are transformable by a simple linear 
equation, 
Cins Per +48 
c 


whereas for magnitudes the equation has the form, 


m’ =m+t+ecC + & 
where e«, wu are "Scale" constants and 4c, 4m are "zero-point" 
constants. For the filter systemvin use for this study, che 
transformation equations are applied to the three indices 
bry, Myr Cy and one magnitude, y which corresponds closest to 


the V of the UBV system. The equations are 
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In their definition of the standard system for "uvby" photo- 
metry, Crawford and Barnes provide a list of carefully studied 
stars with accurately measured values of V, (b-y), M41 and 


c Therefore from an accurate investigation of a group of 
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standard stars, a series of equations can be formulated and 

the unknown coefficients obtained by means of a least squares 
technique. As in the instance of the extinction coefficients, 
mean values for the transformation coefficients are found to 


be sufficient for most observing nights. 


eee DIFFERENTIAL PHOTOMETRY 


The application of both types of corrections is im- 
portant when stars are studied as individual sources by means 
of absolute photometry. However, in the case of variable stars 
certain advantages are realized by applying the techniques of 
differential photometry. Differential photometry relies on a 
comparison of a variable with an immediately adjacent comparison 
star known not to be variable. If the comparison star is close 
@énough in position and colour, then a difference in magnitudes 
as demonstrated by equation 3.lb will, eliminate all first,.order 
exzinckion eG@f£fects.|. sin addition tosthe variable and its com 
Parison eta, periodic» observations of a) check Star are also 
included to ensure that there are no undetected fluctuations in 


the comparison star. The extinction equations become 


MBN-Y 4) = Sede Paces kx for the variable and 
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and if the comparison star has been chosen judiciously so 


that Ax can be-assumed negligible then (b'-y') = (OF ne) Give 
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Simplification of the transformation equations Gan also 
be achieved by the same procedure. If i is the transformed 
value of the visual magnitude for the variable and Si Chae of 


the comparison star then 
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In summary, the use of differential photometry reduces both 


extinction and transformation equations to the following: 
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A complete and detailed paper on the principles of the above 
mentioned corrections together with a variety of methods of 


solution is presented by R. Hardie, 29) 


3,20 TIME 


In the discussion of the eclipsing binary it was 
noted that the two components of the light-curve were the ob- 
served magnitude versus a time scale. Of these two quanti- 
ties, the most accurately measured is time. For the study of 
a variable, the time for each measurement is converted to 
Universal Time, the time at which an observer in Greenwich 
would observe the same event. This 1s accomplished by add- 
ing or subtracting the number of hours ahead of or behind Green- 
wich time as designated by the Standard Time Zone. Considera- 
tion must be given to the implementation of daylight saving time 
where necessary. The hours, minutes, and seconds are then con- 
verted to a decimal fraction of a day. Each day is denoted by 
a running number, specified by the Julian Day, counted from 


January 1, 4713 B.C. The Julian Day begins at noon Greenwich 
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time so that the night hours fall on the same Julian Day. 
(1) LIGHT TIME 


As the Earth moves around the sun in its orbit, the distance 
from the Earth to the variable is constantly changing. For 

a Star on the ecliptic, the starlight reaches the Earth be- 
tween eight minutes earlier or eight minutes later than it 
reaches the sun due to 1ts finite velocity. This cyciie var- 
iation in arrival times has little effect on variables with 
periods of a few days but is readily detectable in short per- 
iod stars. To eliminate this phenomenon a small correction 
is added to the Julian date to convert it to Heliocentric 
time, the time that the light would reach the sun. The mathn- 


(2) 


ematical relation governing the heliocentric correction is 
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A = cos 6 cos a 

Be=Ftan cls wot COS 67 s51in- Oo 
The above relations involve the following quantities 
6, GO, declination and right ascension of the star 

e-angle of the ecliptic 23727" 

R earth - sun distance normalized to one astronom- 

Cal au 
O@ position of earth one thesecliptric Telerred to: the 


vernal equinox 
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The multiplicative constant 0*. 005770 is the time, given as a 


decimal of a day, taken by a beam of light to travel a distance 
OfPone- astronomi calyintt.9) Theinimerical™ values for KR, o7 for 
each day of the year together with the value of the corres 
ponding Julian Day are listed in THE AMERICAN EPHEMERIS AND 


NAUTICAL ALMANAC. 


Soa OBSERVATIONAL TECHNIQUES: 


Prior to any observing, the time as registered by the 
digital clock was set to the standard mean solar time. This 
was first accomplished by the use of a common watch which had 
been checked prior to leaving for the observatory. However, 
with the arrival of the radio, the clock was matched to timing 
Signals broadcast by radio station WWV operated by the United 
States National Bureau of Standards. For work on,EI CEP and 
the comparison star A, the photometer was operated at 1900 
volts. d.c. Bach, filter reading, consisted of. a series. of 
individual voltmeter measurements approximately 12 or more in 
number. A complete set of readings for one star involved 
eight separate filter measurements in the sequence uvbyybvu. 
AiHeers each. titer had. been. positioned, a check was made to see 
if the star under observation was still located in the center 


of the isolating diaphragm; if not steps were taken to recenter 
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it. Therefore an entire series took from three to five min- 
utes to complete depending upon the amount of readjusting of 
the telescope. After each set of data, filter readings were 
taken of the surrounding sky. Included in the observing 


schedule were occasional measurements of the check star B. 


as ee 8 REDUCTION OF OBSERVATIONAL DATA 


Fach set of filter readings was checked to eliminate 
any missed printings or measurements that appeared to have 
been printed incorrectly. The automatic printer recorded two 
voltage readings per second and if two consecutive readings 
for the same time appeared, one reading was eliminated. After 
all necessary corrections had been made, averages were taken 
for each filter reading. Then another average between two 
identical filters were made. From this value was subtracted 
the sky reading for the appropriate filter to give the observed 
intensities. These values for EI CEP and A were then used in 
equation 3.1lb to calculate the observed magnitude change for 
the variable. Since a complete observation of one star took 
a finite length of time, the time of observation was taken as 
the midpoint of the set. Because of the length of the orbital 


period, the preliminary reduction used to plot the light curve 
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CHAPTER 4 


THE OBSERVATIONS 


4.1 The Observations of EI CEP and Comparison Star. 


The observing program for EI CEP which began in Aug- 
ust; 19727 and which was’*términated in*October, 1973; re= 
sulted in 208 individual observations for each of the four 
filters. These represent measurements taken on twenty-five 
separate nights. Although EI CEP was observed on many occa- 
Sions, due to the length of the eclipses and the changing 
nature of the sky conditions, no one minimum was followed 
completely during the course of a single night. A list of 
the differential magnitudes for each filter in the sense 
COMPARISON STAR - VARIABLE together with the differential 
values of the colour index (b-y), and the colour differen- 
and c 


ces m are presented in Table 4.1. The observations 


Vs 1 
are listed in order of the corresponding phase which was 

computed from the values of the period and the epoch of 

Primary conjunction scited oy, Abrami s(TABLE 1.1)... All the 
times of observations were corrected to heliocentric time. 
The estimated values of the errors in the measurements of 
phase and for the magnitudes are listed below the corres- 


ponding title. Graphs of Am versus phase were plotted to 


vyreto the “liqgnt clrve for each I1iter. 
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PHASE 


*.0003 


~0009 
«0020 
©0034 
0042 
0053 
«0078 
-0096 
©0114 
OLLy 
00128 
~0138 
00142 
“ou5> 
00169 
«0183 
0194 
.0198 
O2h7 
0220 
0240 
~0240 
20254 
0263 
00269 
0284 
0295 
20297, 


TABLE 4,1 


OBSERVATIONS OBTAINED ON EI CEPHEI 


U 


5 050 


-1,053 
-1.062 
-1.083 
-0.945 
-1.013 
-0.905 
-0.958 
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(i) The Light Curves 


Upon inspection, the light curves display two well 
defined minima of slightly different depths that occur one 
Nalesan orbital period apart thus indicating Little on ne 
Change in the values of P and Cm calculated by Abrami. Super- 
imposed upon the minima and readily noticeable outside eclipse 
is the presence of a high degree of scatter in the measure- 
ments. This scatter is present in varying degrees in all four 
filters and probably arises from poor sky conditions and poor 
observational technigues. A re-examination of the data and 
the elimination of measurements of two nights of observing 
failed to reduce the scatter appreciably. A night-by-night 
plot of the secondary minima indicated that the scatter was 
to be found randomly distributed among all observing nights. 
No effects due to reflection or ellipticity of the individual 
constituents, noticeable outside of eclipse, could be de- 
tected since the scatter effectively disguised any non- 
eclipse variations. Another feature of the light curves is 
the apparent lower value of Am at the shoulder of the as- 
cending branch of the secondary minimum. In order to ex- 
amine this feature more closely, the individual minima have 


been isolated and expanded and are presented in the next 
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graphs 
(11) The Individual Minima 


The secondary minima are to some extent non-sym 
metrical especially in the y filter. A lack of symmetry is 
also prevelant in the primary minima to a lesser degree 
except for the b filter. The ascending branches of the 
secondary minima appear not to reach the level outside eclipse 
enjoyed prior to the descending branch. It Sneed be 
noted however, that more observations were made around the 
descending branch than the ascending branch. Because of the 
presence of the previously mentioned scatter and a lack of 
additional observations, it is difficult to determine whether 
or not these peculiarities in the minima are physical 
phenomenon caused by presence of gas between the components 
or through some unknown mechanism or are merely a result of 
the tnaccuracy in the data points themselves. [In his 
photo-electric study, Abramz. notes that his light curve 
shows no extra eclipse variation. However this work was 
carried out using an unfiltered photometer and as such is less 
sensitive to any spectral phenomenon that might be associated 
with circumstellar gas. Popper also makes no mention in his 
spectroscopic analysis of any spectral features which might 


be indicative of the presence of any surrounding gas. 
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(111) The Colour Indices 


The next three figures are plots of (b-y) m and 


1 / 
C,against phase. The graph of (b-y) shows the least amount 


of scatter while those of m and c, exhibit progressively 


more. Fart.ot this increase can be attributed to the addition 


Oltmancwersctra factor in both m, and c, which increases the 
targest expected error from, .050Gmag. for (b=y)m to, .100 mag. 


for m, and .110 mag fore Each graph still displays a 


Ll: 
larger degree of scatter than mightebe anticipated 
especially during minima. Values of the mean and standard 
deviation were calculated and listed in TABLE 4.2 for dif 
Ferent sections of the light curve as a test of the amount 
of scatter. For the colour difference Cy the values es the 
standard deviation are somewhat smaller between eclipse, 


whereas for the other two indices the distinction is not as 


easily seen. 
4.2 The Observations of the Comparison Star and Check Star 


As a check on the constancy of the emitted radiation 
from the comparison star (HD 204426, ue = 6.90, GO), measure- 
ments were made periodically on a third star (HD 204408, 


Mv = 6.68, AO). TABLE 4.3 lists the comparison - check star 
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- data obtained during the course of investigation of EI CEP. 
The accompanying graphs portray the magnitude difference 
against the corresponding phase of EI CEP at the time of 

the measurement of the check star. All graphs bear a simi- 
lar degree of scatter as found in the corresponding light 
curves. This fact plus the apparent increase in scatter dur- 
ing eclipse found for the colour indices lends more credence 
to the hypothesis that the distortions in the light curves 
are a result of the observational uncertainty rather than of 


any physical anomalies. 


473.) the Observation®of Standard Stars 

Measurements of the intensity of several standard 
stars were made on September 18, 1972, and October 4, 1972 
from the following list: 

(1) 88 SPEG 

(2) 6 BARI 

(3)% -LOSoTAU 

(4) 8 STAU 

(5) MSie4ORL. 
Upon reduction, it was found however, that the sagging ef- 
fect in the photometer-telescope coupling had rendered the 
measurements unsuitable for accurate analysis. Different 


results were found for observations on opposite sides of the 
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zenith. “Another set of data was taken in 1973 using the two 
stars 6 Lyrae and i Aquilae. From the data obtained on 6 
Lyrae, it was possible to compute the extinction coefficients. 
The other star was observed too low on the horizon to make use 
of the corresponding measurements. The extinction coeffi- 


cients calculated are 


Kae ge Ee 50) 
Ky 20208 §. 0G) 
K, Bee eel Gg?) 
K3 = soe (..055') 


The values in brackets are mean values measured at Kitt Peak 
by Crawford and Barnes. 

Insufficient data made it impossible to calculate 
values of the transformation coefficients and therefore the 
Measurements were left in the natural system during the 
analysis for the orbital elements. Due to the problems 
posed by the observations obtained from EI CEP, it was de- 
cided to maximize the observing time spent on the variable 
and comparison star so no further measurements of standard 


stars were made. 
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CHAPTER 5 


ANALYSIS OF THE LIGHT CURVES AND CONCLUDING REMARKS 


ee THE COMPUTER PROGRAM EVAL 


The basis for the method of Kitamura, discussed in 
Chapter 3, is the evaluation of the incomplete Fourier trans- 
forms from which are computed the characteristic functions 
for that particular system. The evaluation of these trans- 
forms, requiring a numerical integration, was accomplished by 
means of the computer program EVAL written in FORTRAN for use 
with the IBM 36067. The first task required of EVAL was to 
convert the observations from differential magnitudes to lum- 
gnOsicy Units by Jnverting the "equality m ="—255" log &. ~ ihe 


values of the transforms could then be calculated for each 


branch by means of a trapezoidal rule oe of the form: 
1)3he tes ~ _— 
et 8b ‘ : 1 Re : ; * ; 
Sn* = 4 2 E (8, ) Sin n o5 + & (9, ,,) sin nog i+2| oe 
i=1 
Viste le ate estes 
and 
TA 
Cn' = % caG (@,)cos n Beat &' (8,,,)cos n ae) : ee = 


To ensure the best accuracy possible in the integration, the 


program also calculated the corresponding end COurecCwionss: (Viz. 
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Where h is the average interval of the phases of the observed 
points and ¢ is an arbitrary phase angle outside eclipse, 
while %' (0) is the value of the luminosity of the closest point 
near mid-minimum Once the above summations had been carried 
out, it was necessary to obtain a value that would normalize 
the light curve to bring it into line with the theoretical 
values of the light curve obtained from the a -functions which 
were used to calculate theoretical values of the characteristic 
functions. Lhe normalizing factor tor the, lignt curves as 


obtained from the identity; 


K Bs = 5 


En (gy) acy 


The quantities ¢, and 5 represent two different phase angles 
outside eclipse while the denominator incorporates values of the 
incomplete Fourier transforms computed from the observed lum 
tnositkies withn = 0. A normalizing factor was calculated for 
each branch of both minima, with the final value being their 


arithmetic mean. EVAL next computed 
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where A is the normalized value of the depth of minimum esti- 
mated from the light curves and e€ represents the phase angle 
of the point nearest to midminimum. If a value of phase 
angle occurs at exactly mid-ecliopse then EVAL set € automa- 
tically to zero. It is from the above values of LS. and 
L,C. that the program finally calculates the characteristic 


functions’, 


F = iS 

ae 

aa ay ed ee) 
E = LC, 

1- i 


Also evaluated is the ratio ee net which is a slowly 
varying function of the geometrical elements that becomes 
almost invariant in some ranges, and is a convenient com- 
parison against (F,/F >) theo for detecting satisfactory sets 
of geometrical elements. The final task carried out by EVAL 
was to print out tabulated values of the incomplete trans- 


forms, characteristic functions, and the normalized values 


of the luminosity. 
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Sie THE RESULTS OBTAINED BY USE OF EVAL 


Before final reduction, the data were checked 
for measurements corresponding to the same phase. If 
found, the values of the differential magnitudes were 
averaged. This re-examination reduced the number of data 
points to 204. The output from EVAL obtained from the 


observations on EI CEPHEI are listed as follows: 


TABLE (Se. u-filter 
TABLE oe 2 v-filter 
TABLE. Ses b-filter 


TABLE 5.4 y-filter 


The column headed with the letter ,K, is the normalizing 
faerorefor that “particular branen.4 The nextetable tilius= 
trates the values of the estimated depths of minimum with 
their normalized values. The list of the normalized lumin~— 
osities are numerically presented in TABLE 5.6 and graphically 


displaved by eriGe 5.10%5.2, 5.8); 5.4. 


Des CONCLUSIONS CONCERNING EI CEP 

The values of the characteristic functions for 
all four filters display a significant difference between 
branches. A comparison between minima reveals for the 


u.b.y filters that in all instances E < 
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TABLE 5.5 
DEPTH OF MINIMA 


ESTIMATED NORMALIZED 
(MAGNITUDE) (LUMINOSITY) 
PRIMARY SECONDARY PRIMARY SECONDARY 
~1.060 -1.045 649 658 
~0,885 -0,845 639 «663 
~1.065 01,0 sOge 665 


-1.225 -1.155 643 «686 


7 

ft et i 
eeu § FF 
“eee os 


7 
a a 7 : 
; oe 
reves S) a 44 Sc iaihiaal 
: are 
(SeTAee teats raaee evan) 


2 
yandtedis mRNLR ryarosen = YMA 
me ela © M.ta Veter 
coe. of. Pill, O~ ae 
aad, st. no. tf. 


| 6 Ct Sy Bete 


PHASE 


»0009 
-0020 
0034 
0042 
.0053 
-0078 
0096 
0114 
0117 
20128 
0158 
0142 
20155 
.0169 
«OLS 
~0194 
<Ouge 
eO2k7 
0220 
0240 
0254 
.0263 
30269 
0284 
20295 
00297 
0304 
0319 
-0333 
20353 
10559 


TABLE 5.6 


NORMALIZED VALUES OF LUMINOSITY 
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F,sec < F prim, and F.sec < Foprim. The above inequal- 


2 
ities hold for the u-filter if average values between 
branches are taken. This implies)"using the relationships 
given by Kitamura namely foo ae ee and or ees 
that the primary minimum of EI CEP is an occultation while 


Pe secondary 19a transit. This’result 1s sconsistent 


with the conclusions arrived at by Abrami. 


The discrepancy between the characteristic functions 
of the ascending and descending limbs are such as to leave 
Ene spOssipbriity Of a: final solution to the Iqht curves 
somewhat in doubt. While the values do agree to two deci- 
mal places, agreement to the third or fourth is required 
to use the theoretical tables of the characteristic 
functions successfully. This close agreement between values 
is especially required of the ratio CB ere because 
of its slow variation with respect to the geometrical ele- 
ments. The required agreement, however, is not achieved 
as the values obtained for the descending branch of the 
primary minimum are considerably lower than the other three 
for all the filters. Despite the amount of the discrepancy 


between branches, an attempt was made at a solution for 


the b-filter which displayed the closest agreement. An 
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extensive search through the TABLES failed to yield a 
reasonable set of theoretical values to match the obser 


vational values and therefore a complete analysis could not 


be realized. 


5.4 CONCLUDING REMARKS 


This initial test of a narrow passband filter 
system used in conjunction with the University of Alberta 
telescope has posed several problems. First, it appears 
that the employment of this system on a faint star intro 
duces a large scatter. Although the light curve is dis 
Ceotnable: it is not.oLr Sulticient, accuracy: towalilovgton 
complete analysis using the method of Kitamura. A major 
question arises from whether or not the scatter is a re 
sult of simply using too faint a system and overextending 
the bounds of the equipment or rather an intrinsic property 
of a narrow band-small telescope union. It would seem an 
appropriate step to re-examine the uvby system under 
different conditions. This new study should involve an 
eclipsing binary of brighter magnitude, shorter period, and 
Dose byawitie anvexisting: accurate light, curve. The need 
for the increased luminosity is to reduce the back ground 


effects of the sky by increasing the signal from the 
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photometer for all filters yet Maintaining an operating vol-~ 
tage close to the value used on EI CEPHEI (1900V). The 
necessity of a shorter period is to ensure that each eclipse 
may be observed more fully on a single night and that more 
nights of observing are available when the star is in mini- 
mum. The concept of using an eclipsing binary with a prev- 
iously measured light curve is to determine the overall ef- 
fect of the scatter (if any is present) on the shapes of 


the minima. 


If the scatter is still prevalent, then it may be 
necessary to take measurements of the comparison star before 
and, atter the variable to improve the accuracy.) Tis pr0— 
cess should be facilitated with the new mounting which is 
motor driven in both right ascension and declination. Con- 
stant monitoring of the sky conditions is essential so that 
it does not become an appreciable fraction of the voltages 
received from the stars being observed. During observations 
of EI CEP, these readings could be as high as thirty percent 
of the signal obtained in the u-filter. Consideration 
should be given also to the possible extended use of a pulse 
counter to determine if any gain in accuracy can be achieved 


over the digital voltmeter now in operation. 
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The measurement of a bright standard star prior to a 
night's observation may indicate the reproducibility of 

the equipment... If monitored turougnoursthe night, at could 
give some indication of the operation of the photometer 


during the course of the night. 


Although the light curves obtained were inadequate 
for a complete analysis, EI CEPHEI proved an interesting 
system to work with and in itself is worthy of a closer, 


more accurate study. 
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APPENDIX I 


In order to make the most effective use of EVAL 
a short program CORRslipnuld be executed first. The prime 
objectives of this program are to add the heliocentric 
corrections, calculate the colour differences and to order 
the data in a manner compatible with the input requirements 
of EVAL. Since this program together with EVAL was 
written primarily for EI CEP certain changes may be necessary 
to make these programs suitable for other variables. For 
instance to obtain the entire printed and punched output 
from CORR, the programmer will have to delete all the C's 
placed in the first column of some of the statements. These 
were used when the program was rerun and not all the out 
put was essential. The input data for CORR is; 
(1) ‘the titles of the graphs. of the dight curves, 
(b-y),m, , andc, . For a detailed explanation, 
the reader is referred to the COMPUTING CENTRE 
publication entitled SYSTEMS SUBROUTINES in 
particular CGPI: 
C2) the following card contains#the constants: re 
quired by the program; 
(a) RIGHT ASCENSION (ALFA) of the variable 


expressed in radians, F7O, 
(b) DECLINATION (DELTA) of the variable expressed 


in radians, &£7.0 
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ANGLE of the, BCLIPTLG (ECLIP) exexressed 

Bigat cers heli st hon Ue (eC) 

PERIOD of the variable (PERIOD) F9.6, 

EPOCH OF PRIMARY CONJUNCTION (TO) expressed 


by the corresponding Julian date, F13.4. 


next cards put the nightly data in memory. 


furs icard of machinicdht swobsexnvatiens 


contains, 


(a) 


(b) 


(da) 


the) Julman Day ati-the begqinnangrofethe 
observing night (X7DAN) ho. 07 

the number of observations made during that 
Might {NO),i3, 

the normalized EARTH-SUN radius for the 
above Julian Day, (RAD), F9.0, 

the longitude of the EARTH “inj ats orbit 
expressed in degrees (DS), F5.0, minutes 


(SM) > F4.0, seconds, (S5)%,786.0- 


The last two quantities are obtained from tables given in 


the NAUTICAL ALMANAC. Mies lastwcands 24st the weducediaata 


from the digital voltmeter with, 


(a) 


the time of each observation in hours (HR), 


F4.0, and minutes (MIN), F4.0, 
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(b) the differential magnitudes computed from 
the raw data in the order u(xXU), v(XxV), 


b(XB), y (XY), 4F7.0. 


All the above format codes should be considered as flexible 
and easily altered to more suitable values. The values of 
the differential magnitudes for EI CEP were negative but 

for convenience were read in as positive and adjusted within 
CORR.| therefore for data which is) intrinsically positive it 
will be necessary to remove the negative signs present in 
fhe statements that assign these values (XU, KV, XB,;XY) to 


iM TBS lore AV AGHY eeras 31 10) rere lola ga (1) A 


The output data from CORR consists of the printed 


Warues: Of ‘the PHASE, U(l), V(L) > Bs) ally = eth); May 


Cc and I. The value of the total number of data points, 


at 
MAX, together with the PHASE and the four magnitudes are 
punched on cards ready for use with EVAL. The cards are 


punched in increasing order of phase from .0000 to 1.0000 
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with I representing a running number from 1 which corresponds 


to the phase nearest..0000. The. use of I simplifies the 
determination of the limits of the summations used in the 


next program. 
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APPENDIX IT 


phe input catastor EVAL consist of, 

(1) the titles of the graphs of the normalized 
light curves as explained in the preceding 
appendix, 

(2) the estimated values of the depth of primary 
and secondary minima for each of the filters 
expressed in units of magnitude (ULMP, ULMS, 
VLMP, VLMS, BLMP, BLMS, YLMP, YLMS), 8F8.0, 

(3) the following values are the integer numbers, 
I, of the various limits required for the in- 
tegration; 

(a) IFI1, value of phase angle outside eclipse 
(1 form cach, Dranch)s,. la, 

(b) IFI2, a different phase angle further out 
side eclipse (1 for each branch) 14, 

(c) NEAR, the point on the branch being inte- 
grated closest to midminimum (1 for each 
branch), 14, 

(d) IEPS, the point nearest to mid-eclipse on 
any of the two branches forming a minimum 


(1 ftoreeach minimum), 4. 


The values of NEAR and IEPS will be identical for one branch 


in each minimum. The order of the branches being summed 
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over is primary ascending, primary descending, secondary 
ascending, and secondary descending. 

(4) The rest of the deck is comprised of the punched 

Output from CORR. 

It may be necessary at this point to make some corrections 
to EVAL. As in the previous program the magnitudes were 
presumed negative therefore if the data is positive and this 
includes the estimated depth of minimum, it will be neces 
sary to convert these quantities to negative values when 
computing their respective luminosities by changing 
(0,4 * U(IT)) to (-0.4 * U(I) and (0,4 * ULMP) to~(—0.4*ULMP). 
Again as in CORR certain information was suppressed by placing 
a €*¥in column 1, but?’ a quickssurvey of the accompanying 
listing will reveal which statements should be altered. If 
using the program presently stored at the University of 
Alberta, another card containing the anteger 7 ani column 2 
will have to be placed at the end of the input deck. If 
the following listing is being used as a guide no extra 
cards are required. If any unresolvable problems should 
occur, the author of these programs can be reach by address- 
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LORENCE TOMAS ZEWSKI 
Dept. of Astronomy 

Univ. of WESTERN ONTARIO 
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